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Abstract 



Multi-access networks may exhibit severe unfairness in throughput. Recent studies show 
that this unfairness is due to local differences in the neighborhood structure: Nodes with 
less neighbors receive better access. We study the unfairness in saturated linear networks, 
and adapt the multi-access CSMA protocol to remove the unfairness completely, by choosing 
the activation rates of nodes appropriately as a function of the number of neighbors. We 
then investigate the consequences of this choice of activation rates on the network-average 
saturated throughput, and we show that these rates perform well in a non-saturated setting. 

1 Introduction 

Multi-access protocols such as CSMA [IT] as used in the IEEE 802.11 standard have gained much 
popularity for their ability to regulate the access of network nodes to a shared communication 
channel in a simple and fully distributed fashion. A major drawback of these protocols, however, 
is that they can exhibit severe unfairness, in the sense that some of the nodes get starved, while 
others receive good access. One of the main causes of this unfairness is that all nodes are treated 
the same, irrespective of their position in the network. We propose a way of compensating for 
the possible spatial disadvantages of nodes by enhancing the multi-access protocols with local 
information about the immediate neighborhood of nodes. 

Unfairness in wireless networks is an active topic of research. Wang and Kar [IS] considered 
three nodes on a line that only block their direct neighbors, and showed that the middle node is 
starved when the activation rate of all three nodes increases. Such unfairness has been studied for 
more general networks by Durvy et al. [A, 5 and Denteneer et al. [3]. We study the same model 
as in [H 03 |3] , a network with n nodes on a line, in which active nodes block a certain subset 
of other nodes. Unblocked nodes become active, and active nodes deactivate, after exponential 
times. Under these assumptions, the rt-dimensional process that describes the activity of nodes 
is a reversible continuous-time Markov process that possesses an elegant product-form expression 
for its stationary distribution. So this idealized linear network strikes a good balance between 
tractability and complexity, as it retains the essential features of competition and unfairness. 

We assume that an active node blocks the first j3 nodes on both sides, and we say that nodes 
that might block each other are neighbors. Results from [5J [5J [3] suggest that the unfairness 
observed in this model is due entirely to boundary effects. To deal with the unfairness, we shall 
modify the model in one important way: Instead of letting nodes activate at the same rate, we 
allow for node-specific activation rates. By choosing the activation rate Vi of node i as a particular 
function of the number of its neighbors, we can guarantee that all nodes in the network have the 
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same throughput, completely removing the unfairness. Our main contribution is that we prove 
that this fair choice of activation rates vi = v* takes on the extremely simple form 



with "/(i) the number of neighbors of node i, and a any positive constant. We note that these 
node-dependent activation rates are still in line with the distributed nature of the multi-access 
protocol, as v* only requires the number of neighbors, which a node can obtain locally by sensing 
its direct environment. By choosing the activation rates according to |T]) we essentially ensure 
long-term fairness. The first result in this direction is due to Kelly [9], who considered a tree 
network (of which the linear network is a special case) with nearest-neighbor blocking (fi = 1). 

If we choose the activation rates as in ((XJ) and let a —> oo , the throughput approaches 1 / (/3 + 1 ) , 
the highest possible per-node throughput. Recently, Jiang and Walrand 8[, Rajagopalan and 
Shah [14] and Liu et al. [12j proposed protocols for determining the activation rates that can 
attain this maximal throughput in general networks. To this end, they introduced adaptive mech- 
anisms that converge over time to the throughput-optimal activation rates. The main difference 
between our approach and the results in [8l 1141 112j is that we obtain an explicit expression for the 
throughput-optimal and fair activation rates, rather than relying on an algorithm to find these. 
In fact, our explicit solution can be used as a benchmark for validating, in this particular case, 
the various algorithms proposed in [SJ HU [TJ] . Besides giving insight into the operation of the 
network, these explicit rates also provide immediate optimal performance of the network, rather 
than first going through a startup stage during which nodes have to learn the right activation 
rates. The results discussed in our paper are only valid for linear networks, although the Markov 
random field approach from [9] allows us the extend our results to various other topologies, see 
Section El 

The paper is structured as follows. In Section[2]we introduce the linear network in more detail. 
In Section [3] we study some of the key features of the unfairness that arises when all nodes have 
equal activation rates. Section |4] contains the proof of the proposition that the activation rates 
in flj) yield equal throughputs. In Section [5] we investigate the impact of the rates in (p} on the 
network-average throughput, and in Section[5]we discuss the extension to non-saturated networks. 
Section [7] presents some conclusions and further research directions. 

2 Model description 

We consider n nodes on a line, and assume that all nodes are saturated (i.e. have an infinite supply 
of packets available for transmission). Nodes activate according to a /3-hop blocking model, so 
a transmitting node blocks the first j3 nodes on both sides. When node i is blocked, it remains 
silent until all nodes within distance j3 are inactive, at which point it tries to activate after an 
exponentially distributed (backoff) time with mean Node i activates if it is still unblocked 

when the backoff timer runs out. If a node finds itself blocked when the backoff timer expires, it 
waits until all neighboring nodes are inactive once more and then repeats the backoff procedure. 
Without loss of generality, we assume that transmissions last for an exponentially distributed time 
with unit mean. Under these assumptions, the n-dimensional process that describes the activity 
of nodes is a continuous-time Markov process. 

Each state of the Markov process is described as 



where uii — 1 when node i is active. Let Q, C {0, 1}™ be the set of all feasible states. Call ui e £1 
feasible if no two l's in oj are j3 positions or less apart, i.e., w,;Wfc = if 1 < \i — k\ < /3. 

The Markov process that describes the activity of nodes is then fully specified by the state 
space f2 and the transition rates 



(1) 



u = (u)i, ...,oj n )e {o, 1} 
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Here denotes the vector with all zeros except for a 1 at position i. 

Alternatively, we can express the set of feasible states as all states that satisfy a certain system 
of linear equations. Let A be an (n — /3) x n matrix where each row contains /3 + 1 consecutive 
l's, in the following way: 



A = 



( 1 
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Now we can write the state space as fi = {u € {0, 1}" | Aim < C}, where C is the all-1 vector 
of size n. This characterization has a natural interpretation as a set of capacity constraints, and 
nodes can activate only when enough capacity is available. We allocate unit capacity to each 
node, and use the convention that whenever a node is active it uses its own capacity, as well as the 
capacity of all its neighbors to the left. The i-th row of A thus represents the capacity required 
when node i is active. The constraints that arise from the last j3 nodes on the line are redundant, 
and ignoring these leads to the matrix A in 

From the description of the state space as a set of capacity constraints, it is clear that our 
model belongs to the general class of loss networks, see Kelly [10]. Loss networks are known to 
be reversible, and possess product-form solutions. For our Markov process, this product-form 
solution is given by the stationary measure n on D, for which 



7T (w) 







if uj is feasible, 
otherwise, 



(4) 



and where Z n is the normalization constant that makes tt a probability measure. This result is 
well-known in the context of wireless networks, see e.g. [2[ 131 l4[ [T8]. 

Our main concern is with the long-term behavior of nodes, characterized by their throughputs. 
A common throughput-degrading phenomenon in wireless networks is collisions, which may occur 
when multiple nearby nodes transmit simultaneously, causing these transmissions to fail. However, 
we assume that the blocking range (3 is large enough to rule out collisions, so any activity of a node 
contributes to its throughput. Although a model without collisions might seem limited, numerous 
simulation studies show that choosing the blocking range just large enough to preclude collisions 
gives very good performance, see e.g. [3 [TT1 [HI [21]. In fact, the forthcoming paper [16] shows 
that this choice is throughput-optimal when the activation rates are sufficiently large. We study 
the throughput vector 9 = (6>i, . . . , 9 n ), where 9i represents the fraction of time node i is active. 
Denote the total number of feasible states by K, let il = {Sli, . . . , £Ik}, and introduce the n x K 
incidence matrix X such that X^ = 1 when the ith element in the state flk equals 1. Then 



x ■ n, 



(5) 



with n = (7r(Qi), . . . , n(n K j). 

By exploiting the structure of the network, we can construct alternative expressions for the 
throughput in ^ . Specifically, we shall make use of the observation that if node i is active, nodes 
to the left of i behave independently from nodes to the right of i. This leads to the following 
theorem. 



Theorem 1. Define the sequence {Zi)^._ oa such that Zi = 1 for i < 0, and 

Z i = l + v 1 + --- + v h i= l,2,...,/3 + l, 
Zi = Zt-! + ViZi-p-i, i = + 2,j8 + 3, . . . . 

Let the vector of activation rates v — (yx, . . . , u n ) be symmetric. Then 



Zi-B-lZn-i-p 
V% fi . 1 

An. 



1, 



(6) 
(7) 



(8) 
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Proof. By conditioning on whether node i is active, we can decompose the activity of the network 
into two parts, separated by this active node (see [2], Equation (15)), 

e t . = v fi±±^±i^ t (9) 

^l:n 

where Z^j is the normalization constant of a network consisting only of nodes i, . . . , j. For sim- 
plicity we denote Z\ := Z\-,i, and the symmetry of v implies 

Zi-.n = ^l:n-»+l- (10) 

Substituting (fTU)) into (|9|) yields the expression for 0, in ((SJ). By conditioning on the activity of 
node z, we immediately get the recursion relation (0 for the Z;. □ 



3 Unfairness 

We now venture deeper into the issue of unfairness, and assume all nodes to have equal activation 
rates Vi — a . As observed by Durvy et al. [1] and Denteneer et al. [3], the throughput distribution in 
this setting is highly unfair, in the sense that some nodes have a systematically higher throughput 
than others. This unfairness can be explained by the node-in-the-middle phenomenon discussed 
for example in Wang and Kar [TS] and Garetto et al. [8] for the case n = 3, /3 = 1: The middle 
node is in an unfavorable position as it has to wait for both outer nodes to deactivate, whereas 
these boundary nodes only compete for transmission with the middle node. 
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Figure 1: The per- node throughput for j3 = 1 and various values of n and a 

In order to compute the throughput from (|8[). we require information on the Zi. For general 
activation rates vi these cannot be solved explicitly from the recursion ©-([7]) (an exception to 
this are the rates in (fTJ), see Section [3} . However, when all nodes have equal activation rate a we 
can rewrite the Zi as fsee lA")) 

P 

^i = I>iAj-> i = 0,l,..., (11) 

3=0 
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where Xj are the roots of 



\ f3+1 - X 13 - a = 0, (12) 
and the Cj are given by 

A? +1 

(13) 



J (fi + 1)^-0 

Moreover, when the index grows large (i — ¥ oo), the normalization constant (1111) is dominated by 
the largest root A , i.e., 

Z t ~ c X l Q , i -> oo. (14) 

We shall use this asymptotic relation at several places, see (fT6|) and (|56p . 

For ease of presentation, we restrict ourselves to (3 — 1 for the remainder of this section. 
Figures |l(a)f l(d)| show the per-node throughput for various values of n and a. All figures display 



a similar pattern, with the outer nodes having the highest throughput. Moreover, the throughput 
is symmetric in all figures, and exhibits some form of oscillatory behavior. These observations are 
formalized in the following result. 

Proposition 1. For Vi = a > 0, i = l,...,n and j3 — 1, the throughput has the following 
properties: 

(i) Symmetric: 9i = n -i+i, i = 1, 2, . . . , n. 

(ii) Alternating and converging: (— — 9i) is positive and decreasing for i — 1,2, ... , \n/2\ . 
Proposition [1] is proven inlBl 

In Figure Q] we see that for nearest-neighbor blocking, the largest difference in throughput is 
between 6\ and 02- Proposition [Tfii) confirms that this is the most unfair situation, and it persists 
even in larger networks where the node-in-the-middle problem is mitigated by the activity of the 
remaining nodes. In fact, for large networks we have the following result. 

Proposition 2. For Ui = a > 0, i = 1, . . . , n and (3 = 1, 



— , rwoo. (15) 

Proof. We have that d\ cx Z„_2 and 62 oc Z n ^. Using (p~4|) we obtain 

Q 

— — A , 71 y 00 . (16) 

For /3 = 1 we can explicitly solve (fT2l to obtain Ao, and the result follows. □ 
We note that for /3 = 1 , alternative descriptions of Zi exist of the forms 

z^t-^t'+'t^tvq/to), (17) 

where U n {x) is the nth Chebyshev polynomial of the second kind, and 



E 



(18) 



The latter expression can be interpreted as the summation over all possible combinations of nodes 
that can be active simultaneously. 

Figure Q] shows another interesting property of this network. Increasing a leads in many cases 
to a higher throughput for each of the nodes. Hence, in such situations, one may want to increase 
o~ further. However, we also observe that there exists a critical value a* , such that at least one of 
the throughputs Oi decreases as a increases beyond a* . The characterization of this critical value 
is a possible topic for future research. 

Results similar to those presented in this section can be obtained for f3 > 2. As an example, 
Figures 2(a)|2(b) show the per-node throughput for n = 9 and /3 = 2,3, obtained by iterating the 



recursion ([6J)-([7]). Both figures exhibit similar oscillatory behavior as observed for /3 = 1, although 
the oscillation period increases with (3. 
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Figure 2: The per- node throughput for n — 9 and various values of /3 and ct 



4 Fairness 

In this section we present a way to completely remove the unfairness that was addressed in Sec- 
tion [3l In order to do so, we choose node-dependent activation rates vi such that all nodes have 
equal throughput {9\ = 62 = ■ ■ ■ = 9 n ). From (|4]) and ((U we see that in order to meet this objec- 
tive we have to solve a system of n nonlinear equations. It seems that in general this system cannot 
be solved directly. We therefore choose a more indirect approach, for which we first consider two 
special cases that can be solved explicitly. The insight obtained from these exact solutions is then 
used to guess the general solution to the system of non-linear equations. 

The first case is where /3 = n — 2, so that all but the two outer nodes will block the entire 
network. 

Proposition 3. For linear networks with 3 or more nodes, and f3 = n — 2, setting v\ = v n = a 
and V{ — a(l + a) for all other nodes yields equal throughputs 

h= 1 ,, a T\ , i = l,...,n. (19) 
1 + (n — l)a 

Proof. The expression for the throughput in ([S]) can be written as 

01 = Z- X vx{\ + u n ), (20) 
6 i = Z~ l v i , i = 2,3,...,n-l, (21) 
n = Z- l v n (l + vi). (22) 

The inherent symmetry of the model allows us to set v\ =v n . Moreover, for the throughput of the 
other nodes to be equal, we require v% = ■ ■ ■ = v n -i = + V\). If we set v\ — a, and substitute 
this into (|20)) - (|221) . we get a throughput of 



0* = Z- X a{\ + a). (23) 
The normalization constant Z n can be determined by summing over all feasible states: 

n 

Z n = l + y^ j v i + v x v n = 1 + (n - 2)a(l + a) + 2a + a 2 

i=l 

= (l + o)(l + (n-l)a). (24) 
Substituting (gU) into {23]) yields $19§. □ 



The case n = 5, /3 = 3 of Proposition [3] was considered in [3]- The second special case 
corresponds to n = 2((3 + 1), so that a node blocks at least half of the network. 



G 



Proposition 4. For linear networks with n — 2m nodes, m £ IN, and /3 = m — 1, setting 
V{ = a(l + a) 1 ^ 1 for i — 1, . . . , m yields equal throughputs 

i = l,...,n. (25) 



1 + ma 



Proof. To achieve equal throughputs, we see from Q and (J5J) that for the case at hand we should 
solve the system of equations 

v\ + vi(u m+ i H h v n ) = v 2 + v 2 (v m +i H Yv n ) 

= v-i + vz(v m+3 H h v n ) 



= v m + v m v n . (26) 

Indeed, the throughput of node i can be written as a sum over all states in which node i is active. 
Using symmetry, (|26[) can be written as 

V\ + v\(v\ H h v m ) = v 2 + v 2 {v x H h fc'm-i) 

= ^3 + v z {vi -\ h V m - 2 ) 



= v m + v m v 1 . (27) 
Let z^i = a > 0. The solution of (|2"T)l is easily seen to be 

K( = a(l + a) i_1 , i = l,...,m, (28) 

and hence 

9 i = Z- 1 a(l + a) m . (29) 

Summing over all possible states yields 



n m n 



i—l i—1 j—i-\-rn 



m m 



i=l i=l 

= 1 + ((1 + a) m - 1) + mo(l + a) m 

= (l + ma)(l+a) m . (30) 
Substituting ([30]) into j29]) gives ([25]). □ 

It is clear that the complexity of the system of equations governed by ([5]) reduces considerably 
for the choices of j3 discussed in Propositions [3] and |U For general (3 this system remains rather 
complicated, and the direct approach taken in the proofs of Propositions |3] and |4] no longer seems 
to work. 

However, we can use Propositions[3]andH]to make an educated guess about the general solution. 
First observe that the fair activation rates in Propositions [3] and |4] only depend on the number of 
neighbors (nodes within ft hops) that each node has. Denote by the number of neighbors of 
node i, let a > 0, and choose activation rates v* as in (fTJ). We see that this choice is consistent 
with the fair activation rates in Propositions |3] and 2] We now show that v* indeed achieves 
fairness for all ft. To this end, we first show that when the activation rates are chosen according 
to (pj, the recursive relation ©-([7]) for the normalization constant Zi has a closed-form solution. 
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Lemma 1. Let a > and choose v* as in (JTJ) . Then 

Z t = (l + a) 1 , i = l,2,...,n-/3. (31) 

Proof. This can be verified by substituting the solution (|3"Tj) into the relations ([5]) and ([7]). For 

« < P + 1, © gives 

Z, = 1 + a + a(l + a) H h a(l + a) 1-1 = (1 + a) 1 . (32) 

From (0 we see that 

Zi = (1 + a) 2 " 1 + a(l + af{l + a) 1 -' 3 " 1 = (1 + a)\ (33) 
covering the case i > (3 + 2. □ 

With Lemma [TJ we are now in the position to prove our main result. 
Theorem 2. Let a > 0, (3 < n — 1 and choose v* as in |T]). TTien 

Oi = . , " m , *=l,...,n. (34) 
1 + (1 + p)a 

Proof. To prove this result we substitute the normalization constants from Lemma Q] into the 
expression for the throughput in ([8]). We distinguish between different values of i. 
For i > f3 + 1 and i < n — [3 we see that f* = a(l + a)^ and 

VM = (l + ar M , 2 , n _ i _^ = (l+a) B - < -^ (35) 

Similarly for i > j3 + 1 and i > n — /3 + 1 we have i/? = a(l + a) n_I and 

Zt-p-i = (1 + a) 2 -' 3 " 1 , Zn-i.^ - 1. (36) 

For i < (3 and i < n — (3 we have = a(l + a) 1 ^ 1 , and 

= 1, Z n +p = (1 + a)"" 4 -/ 3 , (37) 

and finally for i < (3 and i > n — /3 + lwe have = a(l + a)" - ' 3-1 and 

Zi-p-i — lj Z n -i-p = 1. (38) 

Substituting P5 )l -(|55 jl into © yields 

0i = Z^ofH- a)"-"- 1 . (39) 

We next consider the normalization constant. Choose m such that n = [3 + m. Then 

z n = Z n -i + ^*Z„_^_i 

= Z n -2 + V n _ x Z n -p-2 + V n Z n -p-i 



Z n -P + K+l-i Z n-f)-i- (40) 



Substituting ([3TJ into (00]) yields 



min{m,/3} ^ 



Z n = (l+a) n - 13 + a(l + a) i - 1 (l + a)"-^- i + ^ a(l 

i—1 i=m-\-l 

\n-P-li 



= {l + a) n - p -\l + {(3 + l)a). (41) 
Combining (0TJ» and J3§) leads to ([51). □ 



To better understand why the rates ([T]) successfully ensure strict fairness using only the number 
of neighbors of each node, we study its behavior in the two limiting regimes of light traffic (a 4- 0) 
and heavy traffic (a —> oo). First rewrite ([IJ as 

7(0-7(1) / , 1V . 
3=0 

When a is small, nodes activate slowly, and few nodes will be active simultaneously. In fact, 
the limiting distribution is dominated by states with at most one node active, and node interaction 
(blocking) is negligible. This is reflected in the light-traffic activation rates that follow immediately 
from dHJ: 

v\ = a+( 7 ( i )-7(l))a 2 + 0(a 3 ), a i 0. (43) 

So the first-order light-traffic activation rate v* = a is the same for all nodes. Indeed, if at most one 
node is active (as is the case for a small), there is no blocking, and therefore no need to differentiate 
among nodes. As a increases, states with two active nodes are increasingly likely. This introduces 
node interaction, as nodes may now block their neighbors (all nodes within distance 0). This is 
accounted for in the activation rate by the term (7(1) — 7(I))a 2 , which is linear in the number 
of neighbors. So in light-traffic, only the number of neighbors is of importance, rather than the 
structure of the entire network. This reasoning extends to more general networks. 
When a is large, we see from (|4"2")l that 

v\ =a'rW-7(i)+i + o( Q ,7W-7(i)) j a -> 00. (44) 

Nodes activate almost immediately when they get the chance to do so (i.e. when all neighbors are 
inactive) , and thus the system spends almost all the time in maximal independent sets of active 
nodes. However, as some nodes have a higher activation rate than others, the only states u> that 
have positive limiting probability for a — > 00 are those that maximize the sums of the exponents 
of a over all active nodes, Y^i=A"f(^) ~~ 7(l))I{w 4 =i}) with Iq the indicator function. This is in 
contrast to the heavy-traffic behavior of the system in which all nodes have equal activation rates 
a. Here only states with the most active nodes (maximal independent sets) have positive limiting 
probability when a — > 00. As (fT]) provides fairness even in heavy-traffic, the heavy-traffic fair 
activation rates are such that all nodes are contained an equal number of times in these dominant 
sets. This result strongly depends on the structure of the network, as the maximal independent 
sets may change drastically with the addition or the removal of even a single node. So contrary 
to light-traffic, the simple, locally determined heavy-traffic activation rates (TJ3|) may not easily 
extend beyond linear networks. 



5 Network-average throughput 

The fair rates v* in |T]) are designed to remove the unfairness that arises when all nodes have 
equal activation rates a. In order to compare the two schemes, we want to set their respective 
parameters a and a such that the average per-node throughputs are equal. For given equal rates 
a > we write Zi(a) and 9i(a) for the normalization constant of a network with i nodes, and the 
throughput of node i, respectively. We let 9 n (a) := n~ x Yl7=i ®i{ a ) denote the average per-node 
throughput in a network with n nodes, and we denote 0(a) := limn-j.oo n (a). 

In Section|3]we showed that all nodes have equal throughput a/ (1 + a(0 + 1)) when using the 
fair activation rates. When all nodes have equal activation rates, a closed-form expression for the 
throughput does not seem available. However, we can express the average throughput in terms of 
the normalization constant Z n . 

Proposition 5. Let i>i = a, i = 1, . . . , n. The average per-node throughput is given by 

S n (-) = -^-/z n (a). (45) 
n da I 
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Proof. We have from @ with ui — a that 

n 

8n(cr) = — — - V Zi^p^iZ n -i-p. (46) 

2 — 1 

We compute, using the definition of Zj in Theorem [TJ 

oo n /3 1 2 

^ (^Z^^V^/j)^ = ^(^- r +^GzW) , (47) 

n— 1 z— 1 

where Gz{x) is the generating function of the Z^. Then from (|66[) and some simplifications, we 
get 

oo n 

J2 {Y^Zi-^Zn-i-^ = (1 _ X " (7X , +1)2 - (48) 

n=l t=l ^ ' 

On the other hand, we compute from the explicit form (|66|) of Gz(x) that 

and the result follows. □ 

Using Proposition [S] and the partial fraction expansion (jTTJ) of Z^, we can express the average 
per-node throughput in terms of the roots Aq, • • • \p of (IT^I) . 



Proposition 6. Let Vi = a, i = 1, . . . , n. The average per-node throughput is given by 

a P 
n Q 



(o") = -5, (50) 



where 

P_V A j f n + fi + l \fi + l)Aj \ _ A 3 

^ OS + 1)A,- - ^ V^OS + 1)A,- - jS ((0 + 1)^-0)*)' w f^ w + ^- 

Proof. By d68j) and ([69]) we have 



(51) 



P \ 71+0+1 

z " M =^wrmrTr (52) 

where A., are the (/3 + 1) roots A of (fl"2j) . By implicit differentiation of (fT2|) with respect to a we 
find 

da xfj-ltf + l^-p- W 

Then from ([HJ and JUJ we get 

dZ»(a) = ' / (n + /3 + l)A^ _ (/3 + l)A; + ^ +1 \ dA^ 
da ^ (/3 + i) Aj _ /? (09 + 1)A, - J da 

_^ A " +1 f " + + 1 (1 + ljAj \ 

^ ( / 3 + l)A j -/3V(/3 + l)A J --/3 ((^ + l)A j -^) 2 / 1 ; 

The result follows from substituting ([52)1 and (|54l) into the expression for the average per-node 
throughput (HU). □ 
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When the network grows large (n — > oo) the root of largest modulus, Ao, becomes dominant, 
and tSnj-tSU) simplifies. 



Corollary 1. Let z/j = a, i = 1, . . . , n. The limiting average per-node throughput as n — > oo is 
given by 

e » = W^hrr < 55 > 



Proof. We have, as n — » oo 

(/3 + 1)A -P(p+ 1)A - x ~ ' 08 + 1)A - /? 



\n+l A n+/3+l 

^=7^7T1 o To - TTT\ S (l + °(1))> Q = 7F~T7v\ sC^C 1 ))- (56) 



Hence 



cr P fTA" 3 



(/3 + l)A -p 



and the result follows as crA = Ao — 1 by (fT2j) . □ 

The limiting expression ([55]) occurs in a variety of contexts in [13l [TJ [20l [5j [16] ■ When (3a — > oo, 
the throughput (I55|) simplifies even further. 



Corollary 2. Let ^ = cr, i = 1, ...,n and ief n — > oo. XTie limiting average throughput as 
(3a — > oo satisfies 

((7 ) = ^1_(1 + O (1)). (58) 

Proof. By rewriting (|55l) we have 

^» = ^lI T ^T^-- (59) 

09+l)(A D -l) 

So for 0(a) = 7j£j(l + °(1)) to hold, it is necessary and sufficient that (/3 + 1)(Aq — 1) — > oo. Recall 
from (fT2")l that Ao is such that 

Ag(A -l) = ff. (60) 
Let M > be some positive constant, and assume that /3a < M. Then 

p\%(\ -l)=/3a<M, (61) 

and so (3(Xo — 1) < Af- Conversely, assume that /3(Ao — 1) < -JT for some positive constant -K" > 0. 
Then 

/3a = /3(A - l)Ag < /3(A - l)exp(/3(A - 1)) < Ke K . (62) 

Hence 

/3(Ao — 1) bounded /3cr bounded. (63) 
It follows that a sufficient condition for (/3 + l)(Ao — 1) — >• oo is that (3a — ► oo. □ 

Corollary [2] implies for /3 fixed and a — > oo that (9(cr) — ► g^j- So for n — > oo, both the equal 
and fair activation rates can achieve the maximum throughput by letting a — > oo and a — > oo, 
respectively. This can be explained by the observation that for both sets of activation rates, the 
system spends almost all the time in maximal independent sets of active nodes, thus maximizing 
spatial reuse. 

Next, fix cr > and search for a = a n {a) such that 

n (a) = 



so the network-average throughput is identical for the fair rates and equal rates. For a(a) := 
linin^oo a n (a) we can make this comparison explicit. 
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Proposition 7. We have 



a(a) = A - 1. 



(64) 



Proof. This follows at once by equating (fM|) and (|55|) and solving for a. □ 

It is intuitively clear that imposing fairness may compromise the throughput. From (|34|) it is 
seen that the fair per-node throughputs are bounded above by -g^j, and that this upper bound 
can be approached by letting a — > oo. Corollary Q] shows that, asn-> oo, the average throughputs 
in the fair case and unfair case are equal when a is taken to be Ao — 1. The maximum activation 
rate in this limiting case equals 

ui m = a(l + af = u i , P + l<i<n-f3 (65) 

as is seen from (JIJ. This maximum grows like a since 1 + a = Ao = cr L /^ +1 )(l + o(l)), a — > oo, 
according to (|73l) . Hence, as n — > oo, the fair case achieves the same average throughput with 
activation rates that are of smaller or of the same size as as in the unfair case. The situation 
is considerably more complicated when we keep n bounded. Then it may well occur that 9 n (a) 
exceeds gqrr, which is the upper bound for the throughput achievable by the fair scheme. In 
Figure [3] we have plotted curves a = a n (o~), for the case that n = 10 and various values of j3 in 
the range < a < 20. The occurrence of the asymptotes for the curves with /3 = 2, 5 shows that 
^n(c) > -ppi when a is to the right of the asymptote. Remarkably, in the cases that /? = 1,4,9, 
no such asymptotes exist. This can be explained by the observation that a system operating with 
equal rates a — > oo has at all times T^prl nodes activated, cf. [3]. When is integer, the 
resulting average throughput does not exceed the upper bound -g^j for the throughput in the fair 
case and so no asymptote occurs. 




Figure 3: a n {cr) plotted against a for various values of /3 and n = 10 



6 Non-saturated networks 

The saturated model described in Section [2] assumes all nodes to have an infinite supply of packets 
pending for transmission. We now relax this assumption by introducing queueing dynamics. 
Packets arrive at node 1 according to a Poisson process with rate r, and are routed along nodes 
2, . . . ,n. When node n finishes a transmission, the corresponding packet exits the system. All 
nodes have infinite buffer capacity. These dynamics describe a non-saturated network, where 
packets are forwarded through the network, using nodes 2, . . . , n as relays. As nodes are no longer 
saturated, they may not always have packets available for transmission. When the backoff timer 
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of a node expires, and the node is not blocked, it can activate if its buffer contains at least one 
packet. Otherwise, the node will remain inactive and restart its backoff timer. 

To describe the state of the non-saturated network we have to take into account the queue 
length at each node, as well as node activity. The stochastic process that describes both the 
activity of nodes and their queue lengths is once again a Markov process, but does not have the 
appealing product-form limiting distribution encountered in the saturated network. The added 
complexity of the model is reflected by the lack of available analytical results. When node 1 is 
saturated (i.e. r — > oo) and all nodes have equal activation rate a — ¥ oo the system simplifies, and 
in some cases can be analyzed. For n = 3, (3 = 1 it can be shown (cf. [3]) that both nodes 1 and 2 
are unstable, so that the number of packets in their buffers grows without bounds. Using this, [3] 
shows the throughput of node 3 to be equal to 3/10. Shneer and Van de Ven [T5] demonstrate 
that for general (3 and n — 2/3 + 1, nodes 1, . . . , (3 + 1 are unstable, and the throughput of node n 
is bounded as 1/(2/3 + 2) < n < 1/(2/3 + 1). 

The throughput of node n is of special interest, as it represents the end-to-end throughput of 
the network, that is, the intensity at which packets leave the network. When 9 n = r, all packets 
that enter the network eventually leave. If 8 n < r, packets arrive into the system at a higher rate 
than the network can process, and packets will accumulate at certain bottleneck nodes. The two 
results for the special case in [3] and |15) paint a particularly bleak picture as, despite the high 
activation rate, the end-to-end throughput is well below the theoretical upper bound l/(/3 + 1) 
that can be attained through centralized schedulers. Instead, we investigate the performance of a 
non-saturated network that uses the fair rates ([I]) . 

Figure 2] shows simulation results for the end-to-end throughput in a network using the fair 
rates, with n = 5 and (3 — 1, plotted against the arrival rate r. The dashed line corresponds 
to the network where all nodes have equal activation rate a = 6, and the solid line shows the 
throughput of a network with fair activation rates {!]) and a = 11.68. The values of a and a 
are chosen such that the average per-node throughput is equal for both activation schemes, see 
Section [5j The network with equal activation rates performs poorly. When the arrival rate grows 
beyond a certain threshold, node 2 saturates and the throughput drops. The network with fair 
rates, on the other hand, can sustain higher arrival rates and experiences no throughput drop 
when in overload. In fact, the end-to-end throughput approaches the per-node throughput in 
the corresponding saturated network (indicated by the dotted horizontal line) . So the network is 
stable whenever r < a/(l + a(f3 + 1)). 



85 




r 



0.2 0.4 0.6 0.8 



Figure 4: The end-to-end throughput of a network with equal activation rates (dashed) and fair 
activation rates (solid), plotted against the arrival rate at node 1 

We see from Figure [4] that as r increases, the end-to-end throughput eventually reaches the 
per-node throughput in the saturated network for a = 11.68. Sustaining this high end-to-end 
throughput requires an efficient activation scheme, in which maximal independent sets of active 
nodes are persistant for a long time, yet all nodes get equal access. Figure [5] shows a typical 
sample of the activity pattern, for both the equal rates and the fair rates. The arrival rate equals 
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r equals 0.47, so the network with fair rates is stable. The barcode plots indicate for each point 
in time whether nodes are active (black) or inactive (white). We see from Figure 5(a) that for 



equal activation rates, the pattern is very irregular. Node 1 is more active than the other nodes, 
indicating that the network is in fact unstable, and node 2 is the bottleneck node. The fair 
activation rates result in a more regular pattern, and maintain certain maximal independent sets 



of active nodes for a long period of time. Figure 5(b) shows that the maximal independent sets 



{1, 3, 5} and {2, 4} are persistent. So instead of sending a single packet before relinquishing access 
to the medium, nodes send packets in large batches. 
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(a) Equal activation rates 



(b) Fair activation rates 



Figure 5: The activity charts for both equal activation rates and fair activation rates for n = 5, 
(3=1 and r = 0.47 

The length of the dominant periods of a maximal independent set of active nodes, and therefore 
the size of the batches, is related to the queue lengths of the nodes in this set. Typically, when one 
of the queues of the dominant maximal independent set empties, dominance reverts to the other 
maximal independent set. Figure [S] shows the queue lengths of all nodes over time. It shows large 
fluctuations in queue length, which are directly related to the large batch sizes. Nodes 4 and 5 in 
particular empty relatively often, suggesting that these are in fact the nodes that typically break 
the dominance of their own maximal independent set. 
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Figure 6: The queue lengths for the various nodes with fair activation rates, n = 5, j3 = 1 and 
r = 0.47, on a large time scale 

In contrast to the saturated network, the non-saturated network may waste capacity when 
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unblocked nodes cannot activate because of empty buffers. Figure [7] is similar to Figure [HI except 
for the smaller time scale. It suggests that queues rarely empty, and moreover do not stay empty 
for a long time, so the wasted capacity is very small, if not negligible. 
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Figure 7: The queue lengths for the various nodes with fair activation rates, n 
r = 0.47, on a small time scale 



5. 



1 and 



The results presented in this section carry over to general n and f3. The system with fair activa- 
tion rates can sustain arrival rates up to the per-node throughput of the corresponding saturated 
network, and maintains its throughput when r grows even larger. The network with equal activa- 
tion rates with the same (saturated) average per-node throughput fails to attain this end-to-end 
throughput. Whenever this system is overloaded, its throughput drops. Recall from (|34| that 
for a — > oo, the saturated network approaches a per-node throughput of 1/(0 + 1). So the cor- 
responding non-saturated network even approaches the theoretic upper bound on the end-to-end 
throughput. 

The high end-to-end throughput of a non-saturated network using the fair rates, and the poor 
performance of such a network where all nodes have equal activation rates, can be explained by 
the result derived for saturated networks in Sections [3] and |H Indeed, although non-saturated 
networks are more complex than their saturated counterparts, the basic premise that nodes in the 
middle of the network have an unfavorable position compared with nodes on the network boundary 



7 Conclusions and outlook 

In this paper we studied the unfairness in linear multi-access networks. We proposed node-specific 
activation rates in as a function of the number of neighbors, and showed that these rates provide 
equal throughput for all nodes. The rates increase with the number of neighbors. Intuitively, this 
structure can be explained by the observation that, as the number of neighbors increases, a node 
needs a higher activation rate to retain its throughput. Consequently, the rule in (TTJ, which is 
exact in linear networks, might serve as a heuristic in more complex networks. 

Finding exact expressions for the activation rates that provide strict fairness for networks 
beyond the linear network is challenging. Kelly [9] obtained results for trees with nearest-neighbor 
blocking, and finds that rates such as in (JT]), where nodes on the leafs of the tree have lower rates 
than those in the stem of the tree, provide strict fairness. For such trees, it seems possible to 
extend Kelly's result to the /3-hop blocking situation. Also other regular topologies such as certain 
grids appear to admit similar analysis. 
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Finally, let us mention that regular networks like grids or trees may not always be a good 
representation of topologies encountered in practice, which in general are less structured. The 
results obtained in this paper, however, rely heavily on the diagonal structure of the capacity 
matrix A in @, which only exists for certain well-structured networks. For more general networks, 
and hence more general matrices A, the objective of equal througputs boils down to solving the 
system of nonlinear equations that follows from ([5]). In fact, ([5]) can be described in terms of the 
the mapping (with v — {y\, . . . , v n ) £ (0, oo) n ) 

V ^e(v) = z n -x-u = ( J2 fie) 

V uen i=i J k=i n 

with 9{v) £ (0,oo)". It can be shown that the mapping 6 is globally invertible on (0, oo) n . Thus, 
given a vector c <E (0, oo) n , there is a unique v — v{c) £ (0, 00)™ such that 6{v) = c. When c has 
identical entries, this corresponds to all nodes having equal throughput. The full analysis of this 
fixed-point equation is rather involved and will appear elsewhere. 



A Normalization constant 

We take a closer look at the normalization constant Zi defined in ©-Qj in the regime where all 
nodes have equal activation rates Vi = a. This section closely follows Pinksy and Yemini [13] and 
Janssen et al. [T5J. All proofs can be found in [TB] . 

The generating function Gz (x) of the Zi is given by 

00 1 1 8+1 

G z {x) = ^Z lX = {x _ m _ x _ ax3+1 y (66) 

Let Ao, ■ • • , Xp denote the (3 + 1 roots of 

\' 3+1 - A' 3 - a = 0. (67) 

These roots can be shown to be distinct, and from Rouche's theorem it follows that there exists 
a unique positive real root Ao such that Ao > |Aj|, j = 1, . . . , f3. We can obtain the Zi from the 
generating function by applying partial fraction expansion, which gives (|16j. Proposition 1) 

a 

Z i = s £ j c j \\ ,i = 0,l,..., (68) 
j=o 

where \j are the roots of (|67|) . and 

When the index grows large (i — > 00), the normalization constant is dominated by the largest 
root Ao, thus simplifying considerably, i.e., 

Zi = coA (1 + o(l)) , i ->oo. (70) 

The roots Ao, . . . , A^ have representations in terms of infinite scries expressions. Let (x) n — 
T(x+n)/T(x) denote the Pochhammer symbol. Depending on the value of cr, either of the following 
two expansions applies (See [TB], Proposition 7 and 8). For small o > 0, we have 

A (cr) = l + 2^ jj ® ' ( 71 ) 

1=1 

A,-W = E^§^^, 3 = 1,2,...,/?, (72) 
1=1 
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where Wj = a 1 ^ 4 ^ 1 ^ x / 2 )/P . The series expansions in (I7T1) and (I7"2"j) converge for 



08 + 

and diverge otherwise. For large <r > 0, we have 



-I 
0+1 



A^^E^TT^ (73) 



1=1 



where Vj — (T 1 /(^+ 1 ) e 2 '™-' 'I '(/ 3 + 1 ). The series expansion in (|73|) converges for cr > and diverges 
otherwise. 

B Proof of Proposition 1 

We first establish an auxiliary result. Define a(i, l,n) as the number of states in which exactly Z 
nodes are active, including node i. For successive nodes, the following relation holds. 

Lemma 2. Let n € IN, i < [§] - 1. Then: 

a(i,l,n) = a(i + 1,1, n), I < i, (74) 
a(i,l,n) > a(i + 1,1, n), i odd, i < I < \n/2] , (75) 
a{i, I, n) < a(i + 1, 1, n), i even, i < I < [n/2] . (76) 

Proof. The proof is by induction on i. Conditioning on activity of node 1 and node n yields the 
relations 

a(i,l,n) = a(i — 2,1 — 1, n — 2) + a(i — l,l,n— 1), (77) 
a(i, I, n) = a[i, I — 1, n — 2) + a(i, l,n — 1), (78) 



with boundary conditions 



a(0, 1, n) — for all n and I; 
a(l, I, n) = 1 for I > and all n; 
a(l, I, n) = for I < and all n. 

Hence, the initialization step of the induction is 

a(0, l,ri) < a(l,l, n), < I < [n/2] , 
a(0,l,n) = a{l,l,n), l<0. 

Consider odd i < [n/2] — 2, let i + 1 < I < [n/2] , and assume a(i, I, n) > a(i + 1, 1, n). Using ([77)) 
and |J7J]) we get 

a(i + 1, 1, n) = a(i + 1, 1 — 1, n — 2) + a(i + 1, Z, n — 1) 

< a(i,l- l,n-2) + a(i + 1,1, n- 1) = a(z + 2,Z,n). 

This proves assertion (|75p . Assertions (|74l) and (I7B1 can be proved in a similar manner. □ 

We now use Lemma [2] to prove Proposition [T] 
Proof. (Proposition [1} Assertion (i) can be shown by rewriting the throughput as follows: 

6i = Z- 1 E a(i, I, n)a l = Z' 1 ^ a ( n - * + 1, ') n W = ^n-i+i- (79) 
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To prove assertion(ii) we first show that (— — 6j) is positive. That is, 

(-iy(6 l+1 - 9i) = (-lyZ- 1 Y, (<*(* + 1, h n) - a(i, I, n)) a 1 

i 

Ln/2j 

= (-lyz- 1 (°(* + z > n ) - a (*> z > n )) ° l > °> ( 8 °) 

i=i+l 

where the inequality follows from Lemma [5] Using (|80p . Proposition [TJii) follows from 

{-i)\e l+1 - 9t) = {-iy (e i+1 -z- 1 ^ <*(», i, ™y ) 

= (-1) 1 - Z- 1 5^(o(i, I - 1, n - 2) + o(i, I, n - 1))ct ; ) 
> (-1) 4 - Z" 1 ^(a(i, Z - 1, n - 2) + a(i + 1, /, n - l))cr*) 

- (-ir (^+i - z- 1 E < i + 2 > *> n V) 
= (-i) 4+i (0 4+2 -^ +1 ). 

This completes the proof. □ 
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